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Abstract 

The light-front dynamics (LFD) of the scalar field model theory is analyzed to solve the 
two-body bound-state problem. The light-front two-body bound-state equation is extended to 
the full LFD kernel including the ladder, cross-ladder, stretched-box, and particle-antiparticle 
creation/annihilation effects to study the contributions of higher Fock-states. The light-front 
two-body equation is also modified by the term corresponding to the self-energy corrections and 
counter-terms. Using the variational principle, we obtain the numerical result of the binding 
energy B versus the coupling constant a for various mass ratios of the constituent particles 
including the cases of non-zero exchange particle mass. We also discuss the correlation between 
the mass spectrum and the corresponding bound-state wavefunction. 

1 Introduction 

For an accurate calculation of the spectra and wavefunctions of hadrons, it is important to include 
the fundamental relativistic effects such as the correct relativistic energy-momentum relation, re- 
tardation effects, and particle-antiparticle creation and annihilation. Since the quarks and gluons 
inside the hadrons have negligible masses and interact very strongly among themselves via quantum 
chromodynamics (QCD), the relativistic effects play a significant role in analyzing the mass spectra 
and the wavefunction related observables (e.g. form factors, generalized parton distributions, etc.), 
in particular, for the low-lying hadrons. Although the ultimate goal is to analyze the bound-state 
problem in QCD, prior to getting into the real complicated nature of QCD, we may first investigate 
much simpler bound-state system provided by the Wick-Cutkosky model [l|,[2| and the similar scalar 
field model theories. 

C. Savkli et al. [3j-[5j have already used a powerful numerical approach known as the Feynman- 
Schwinger representation approach and investigated the scalar field model theory which they called 
X 2 </>-theory within the given precision of the numerical computation. The earlier work with the 
Feynman-Schwinger representation to scalar-scalar bound states was presented in Ref. [6| . The works 
of C. Savkli ct al. 3-5] investigated not only the stability of x 2 0-theories but also the cancellation 
among the vertex corrections, overlapping self-energy, vacuum polarizations, etc. In order to get 
the full result for two-body bound-states up to the second order in the coupling constant a, it has 
been discussed that it may already be a good approximation just to include the ladder, cross-ladder, 
stretched-box, and the relevant self-energy corrections and counter-terms. This encouraged us to 
look into the scalar field model theories with the analytical tools available to us. 

An analytically tractable tool for the bound-state problem and also known as the most orthodox 
tool for dealing with the relativistic two-body problem in quantum field theory is the Bethe-Salpeter 
formalism utilizing the Green's functions of covariant perturbation theory. In the ladder approx- 
imation of the Bethe-Salpeter formalism the bound-state [l], 0] problem has been analyzed for the 
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system of two particles interacting with a third scalar particle. However, this approach has funda- 
mental difficulties with the relative time dependence and in systematically including higher-order 
irreducible kernels such as crossed diagrams and vacuum fluctuations [8]. 

An alternative approach which can remove these difficulties and restore a systematic perturbative 
calculation for obtaining higher accuracy is the reformulation of the covariant Bethe-Salpeter equa- 
tion at equal light-front time, r = t + z/c |9|-ll if . This is equivalent to expressing the Bethe-Salpeter 



equation in the infinite momentum frame [12h17] . The light-front quantization method [181 . |19j pro- 
vides a relativistic Hamiltonian formalism and a Fock-state representation at equal light-front time 
t, which retains a lot of the simplicity and utility of the Schrodinger nonrelativistic many-body 
theory [18| . This method not only suppresses the vacuum fluctuations but also systematically in- 
cludes cross diagrams when higher Fock-state contributions are taken into account. The relativistic 
Hamiltonian dynamics with equal light-front time r has been known as the light-front dynamics 
(LFD). 

Relativistic two-body bound-states have been anal yze d with the light-front formalism of the 
Bethe-Salpeter approach in the Wick-Cutkosky model 2(J 21 1. The light-front ladder approxima- 



tion in the Wick-Cutkosky model has been extended to the lowest order light-front Tamm-Dancoff 
approximation, which includes the self-energy corrections and counter-terms [22T |. The cross-ladder 
and stretched-box up to second order in the coupling constant a have also been included by V.A. Kar- 
manov et al. [23|, [2J|. However, the higher Fock-state contributions due to the particle-antiparticle 
creation/annihilation process have not yet been included in the cross-ladder contributions |25l. l26j|. 
Also, the light-front bound-state analyses in the scalar field model theories have largely been limited 
to the Wick-Cutkosky model which describes the bound-state system with the equal constituent 
mass and the zero mass of the exchange particle. Therefore, we investigate the contributions of 
the higher Fock-states due to the particle-antiparticle creation/annihilation process and analyze the 
two-body bound-state problem in various combinations of masses for the constituent and exchange 
particles: e.g. mi ^ mi and A ^ 0, when the two constituent particles <f>\ and <p2 have the masses 
to i and 777,2, respectively, and the exchange particle has the mass A. 

In Section 2, we show the light-front formalism of the Bethe-Salpeter equation with the variational 
method and provide our variational wavefunction suitable for the case of non-zero exchange particle 
mass A. In Section 3, we present the numerical results of the spectrum calculation and relate them 
to the wavefunction renormalization and the higher Fock-state contribution. The conclusion follows 
in Section 4. 



2 Formalism 

2.1 The Bound-State equation in LFD 

For simplicity, we consider the scalar field model theories which describe bound states of two scalar 
particles <j)i,4>2 with masses mi and TO2 exchanging another scalar particle x with mass A. While 
mi =/= 7772 and A ^ in general, the model with mi = mi and A = has been known as the 
Wick-Cutkosky model [1, 2j. The relevant interaction lagrangian is given by 

£ = fl(&& + (2-1) 

where g is the coupling constant with the dimension of mass and (f>i (i — 1,2) is the conjugate field 
of tpi. The Bethe-Salpeter equation of this type of model has been reformulated in LFD and the 
light-front ladder approximation has been extended to the lowest order light-front Tamm-Dancoff 



approximation including the self-energy effect 22J . The cross-ladder and stretched-box up to second 



order in the coupling constant a have also been included [23l. \2A\. 

In this work, we also go beyond the light-front ladder approximation to include the cross-ladder, 
stretched-box and the higher Fock-state contributions. In particular, the higher Fock-state con- 
tributions (the particle-antiparticle creation/annihilation effect), which has not been numerically 
analyzed in the previous literature [24T]26j |. are incorporated. The light-front bound-state equation 
for this theory is given by 
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x 1 — x 

dy d 2 lj_ 



y(l - y) 16n' c 



16tt 2 

K(x, k±;y,l±)%l)(y,l_ L ) , 



f(x,k±) \ ip(x,kx) 



(2.2) 



where K(x, k±] y, l±) is the kernel of the bound-state equation, f(x, k±) is the self-energy correction, 
tp(x, k±) is the light-front wavefunction of the bound state, and M is the mass of the bound state. We 
denote the light-front longitudinal momentum fraction as x, y, etc. and the transverse momentum 
as k±, l±, etc. The kernel K(x,k±;y,l±) is provided up to order g 4 including the ladder (L), 
stretched-box (SB), cross-ladder (CL), and higher-Fock (HF) kernels: i.e. 

K{x, k±;y, F_l) = g 2 V L (x, k±;y, T±) 



dz 



16tt 3 



V SB (x, k x ;y, l'x; z,jx) + V UL (x, k ± -,y, l'x;z,jx) + V Ht \x, k x ;y, H;z,j±) 



rCL 



HF, 



(2.3) 



while the self-energy corrections including counter terms [22j are given by 



f(x,k x ) = - I dzlog 



\ 2 z + m 2 (l- z) 2 



x (1 — x) 
mi m 2 



(2.4) 



As the works of C. Savkli et al. indicated the cancellation among the vertex corrections, 

overlapping self-energy, vacuum polarizations, etc., we do not include them in the present work 
but limit here only up to the four-body Fock states in LFD. We will be content with just taking 
into account the kernel K(x,kx',y,lx) in Eq. (|2.3j) and the self-energy correction term f(x,kx) in 
Eq. (|2.4j) . While the ladder kernel (V L ) is given by 



V L {x,kx,y,lx) 



x <-> y 
k x «-> Zj 



{x - y ) 



(x -y) M 2 _ jj+m! _ (k±-T±) 2 +X 1 _ kl+ml 



(2.5) 



the stretched-box (V SB ), cross-ladder (V CL ), and higher-Fock (V^) kernels can be written in the 
form: 



V SB (x,k x ;y,lx;zJx) = 

V CL (x,kx;y,lx;z,jx) = 
V HF {x,k ± ;y,lx;z,jx) = 



z(x -z)(z- y)(l - z) Vl 
1 

z(x — z){z — y)(l — x — y + z) 
1 

z(x — z){z — y)(l — x — y + z) 



8=3 

[V A + V B ] 



(2.6) 

(2.7) 
(2.8) 



where V\, V2, V3, V4, V5, Vq, V7, V%, Va, and Vb are presented in Appendix A. The corresponding dia- 
grams are shown in Figs. |T|— [4] In particular, the particle-antiparticle creation/annihilation process 
which was not included in the previous cross-ladder analysis [24{ is shown in Fig. |4] We denote this 
contribution as V HF . 
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Figure 1: Ladder LFD graphs. 
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Figure 2: Stretched Box LFD graphs. 



2.2 Spectrum Calculation with Variational Method 



In order to solve Eq. (12.2[) numerically, we utilize a variational principle for a dimensionless cou- 
pling constant a given by a — g 2 /(lGirm ^2). Taking the expectation values with a variational 
wavefunction in Eq. (|2.2j) . we get the quadratic equation in terms of a: 



M 



kl 



l-x 



a (y L ) + aip^f) + a 2 (V SB + V CL + V HF ) , (2.9) 



where more explicitly each of the expectation values are given by 

dx d 2 k± +/ 

in I / 



, k^+m 2 k 2 L + m\ 



M z - 



l-x 



x(l — X) 167T 3 



^(x,k ± ) i M 



k 2 , + m? k 2 + m 



1 -x 



1p(x,kx)- 

(2.10) 



dx d 2 k 

x(l — x) 167T 3 



— tp^(x, k±) — - — -f(x, kj_)ip(x, kj_) . 



IT 



(2-11) 



dx d 2 k± dy d 2 l± 
y) 167T 3 



(V L ) = 16^m 1 m 2 / uy -^ ^(x,k ± )V L (x,k r ,y,l±MyJx), (2.12) 

J x{l — x) lD7r J ?/(l ' ' — 



(V ) = (16 7 rm 1 m 2 ) J ^^J^ J ^y)l^ J dz ^^ 

x V SB (x,k±;y,l±;zJ±)ip(yJ±), 



(2-13) 
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Figure 3: Cross Ladder LFD graphs. 




Figure 4: Higher Fock LFD graphs. 
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(V HF ) = {l&itrmmzY 



dx 



dy d 2 lj 



dz ■ lj^{x, fcj 



x(l — x) 167T 3 J y(l — y) 167T 3 J 167T 3 

x V HF (x,k±;y,l±;zJ±)i/j(y,l±) 



(2.15) 



For fixed binding energy B = mi + mi — M, we solve the quadratic equation of the coupling constant 
a and minimize the corresponding expectation values to determine the optimum relation between 
the binding energy B and the coupling constant a: 



k^+mi 
1-x 



2 (yCL + ySB + 

(2.16) 

The minimum of a is found by varying the parameters in the trial wavef unction if)(x,k±_). A judicious 
choice of trial wavefunction is important to get the closest result to the true minimum of a. To 
achieve this goal, we consider available exact solutions in some limiting cases, e.g. the solution of 
the non-relativistic Schrodinger equation for the Coulomb interaction. We also take into account the 
relation of the light-front bound-state equation to the covariant Bethe-Salpeter equation, since some 
analytic spectral functions for the solution of the Bethe-Salpeter equation are known in both weak 
and strong binding limits [27j |. For the 15 state, we take the variational wavefunction parameterized 
I'.v 



ip(x, k_i 



C 2 



x X—x 



c 2 



k 2 ± +m'{ k 2 ± +r 



-4A A+ .h(ra\ + mf) 



(ZA2 



(m? — m|) 2 
3 P - ^ 



c 2 



(l + |2x-l| 
(2.17) 



where the normalization constant A^i s cancels in the ratio given by Eq. (|2.16[) and C is the variational 
parameter. The factor (1 + \2x — 1|) in the denominator of Eq. (12.171) stems from the weak binding 
spectral function of the Bethe-Salpeter solution as shown in [28]. Eq. (|2.17p is the variational 
wavefunction suitable for the case of non-zero exchange particle mass A. We tried the variational 
wavefunction without the A term in Eq. (|2.17j) as well as other trial wavefunctions different from 
Eq. (|2.17j) . but found that the obtained expectation values of a were all larger than what we get 
using Eq. (|2.17[) . This convinced us that the wavefunction given by Eq. (|2.17p is an improvement 
over any other variational wavefunctions that we have considered 



may 



31] 



The origin of Eq. (12371 

be attributed to the exact solution of the Schrodinger equation with the Hulthen potential [2 
Although the Schrodinger equation with the Yukawa potential has not yet been analytically solved, 
the Hulthen potential behaves like the Yukawa potential for small values of r and the exact solutions 
have been obtained for the Schrodinger equation with the Hulthen potential. With the form of the 
variational wavefunction given by Eq. (|2.17l) . we are well equipped to compute mass spectra even 
for the case of A ^ 0. 



2.3 Analytic Calculation of Spectrum and Non-Relativistic Limit 

A simple analytic relation between the coupling constant a and the binding energy B may be 
attained by following the method presented in Refs. 2o|, 22 1. For example, if we take C = M in 
Eq. (|2.17p for the case of mi = m-2 = m and A = 0, we get 



ip(x, fcjj = — 



N 



(2.18) 



M 2 



x(l- 



(l + |2z-l|) 



G 



Without the factor (1 + \2x — 1|) in the denominator of Eq. (|2.18p . this wavefunction corresponds to 
the IS state wavefunction of the Schrodinger equation with the Coulomb potential. Eq. (|2.18[) was 
used to find the analytic expression for the binding energy of the IS state of the ladder approximation 
originally in [2(J [2l| . Including the self-energy correction, the following formula was obtained [22J : 



- = {(«- 2)/(« - l) 1/2 }{^ - taa- 1 ^ - l)- 1 / 2 )} + log(4/u) + ^ log 1 , 



u-2 



(2.19) 



where a = g 2 /(167rro 2 ) and u — m 2 /(m 2 — M 2 /4). In the zero-binding limit, u — > 00, Eq. (|2,19[) 
leads to — = It-v/m, which is the well-known Balmer formula for the IS state: i.e. 

a 2 v 5 



2m- M = ma 2 /A . 



(2.20) 



Similarly, following the method presented in |22| as well as in [20j, we obtain the following relation 
for mi = 77i2 = m and A ^ 0: 



= 



dt 



(i-2)(l + A^) 2 l gj 




u[4+t(t-4)^-] 




4(t- 


-4)+u 


4+(t-2)(t-4)^£- 


r) 


(&u-2) [t(t - 4) +4u] 


t(t 


-4)(1 + A^) 2 + 


4// 



dz log 



2z(l - z) 



A 2 



(2.21) 



2~Z + u(l — z) 5 



For A = 0, Eq. (12.21)) can be integrated and reduced to Eq. (|2.19[) . We will discuss these relations 
further in the zero-binding energy region when we present our numerical results in Section 3. 



2.4 Wavefunction Renormalization and Probabilities of Lowest and Higher 
Fock-States 

The light-front quantization method provides a description of a bound-state \B) in a Fock-state 
representation at equal r: 



\B) = + {H>x\B)\H>x) + {WXX\B)\WXX) 



(2.22) 



where the light-front two-body wavefunction ip(x,k±) corresponds to the lowest Fock-state ampli- 
tude (4>cf>\B). The probability to find the two-body state \</)<j>) is given by the integration of the 
wavefunction squared: 

dx d 2 k± 



x(l — x) 167T 



-\^{x,kA_)Y 



(2.23) 



where the subscript "Low" of Plow is introduced because the two-body state is the lowest Fock- 
state. The probability to find the higher Fock-states can be obtained by differentiating the kernel 
K(x,kx;y,lx) m Eq. (|2.3[) with respect to M 2 as discussed in Refs. [22, 32 1 : i.e. 



dx 



kj_ dy d 2 l 



c(l — x) I671" 3 y(l — y) 167T 3 



£,/,t 



tjy[x, k±) 



f dK 
\ dM 2 



(2.24) 



where the subscript "High" of Pmgh indicates the contributions from the higher Fock-states such as 
three and four bodies through V L , V s B , V CL , and V HF . Including the self-energy correction, we 
may also renormalize the two-body wavefunction as discussed in Ref. (28l |: i.e. 



dx d 2 k± 

x(l — x) 167T 3 



,am\mi df(x, k±) 
' n ' dM 2 



\^{x,k ± f 



(2.25) 
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where ip(x, k±) is the wavefunction including the self-energy correction. The probability to find the 
higher Fock-states including the self-energy correction is given by 



Pm g h 



dK \ 
dM 2 / 



dx d 2 k± dy d 2 l 



x(l — x) 167r 3 y(l — y) 16ir 3 



4>'{x, k±) 



f dK 
[dM 2 



(2.26) 



Since the kernel K can be decomposed into the contributions from V L , V s B , V CL , and V HF , the 
probabilities Plow, Plow, Pmgh, and Pm g h can also be decomposed into the corresponding contribu- 
tions. For example, up to the ladder and cross-ladder, we may define the following probabilities): 



P, 



P L = 

Low 



L+CL 
Low 



Pi 



pL+CL 
Low 



pL 
-'High 



dx 


d 2 k ± 


x(l - 


x) 16-7T 3 


dx 


d 2 k ± 


x(l - 


X) 167T 3 


dx 


d 2 k ± 


x(l - 


x) 167T 3 


dx 


d 2 k ± 



\^L(x,k ± )\\ 

\iPl+cl(x, fc±)| 2 , 

amim 2 . df(x, k±) 
7T ' dM 2 

amim2 . df(x, k±) 



P 



L+CL 
High 



pL 
-'High 



X(l — X) 167T 3 

d(aV L ) 
dM 2 

dx d 2 k± 

x(l — X) 167T 3 

d(aV L +a 2 V CL ) 
dM 2 

dx d 2 k± 

x(l — X) 167T 3 

dx d 2 k± 

x(l — X) 167T 3 

d(aV L ] 
dM 2 



dM 2 



dy d 2 l 



\^l{x, fcj_)| 2 , 
\^l+cl{x, k±)\ 2 , 

dV L 



^4(a!,Ax)(^)^(»,r ± ), 



2/(1 - j/) 16 



^ d2 ^ ^ ( x t. 
2/(l~y) 167^+^'^ 



( dV L 



\ dM 2 



fJ 



dy d?lj 



d 2 J± 



fdV CL \ 



= -.9 



P 



L+CL 



High 



dx d 2 k± 

x(l — x) 167T 3 

d(aV L + a 2 y CL ) 
dx d 2 k± 

x(l — X) 167T 3 

dx d 2 k± 

x(l — X) 167T 3 



c?2/ d 2 Z 



2/(1 - y) 16tt 3 



-^(x,fci) 



(dM 2 



dy d 2 l 



167T 3 ^ L+CLV 



<9V L \ - 



V'l+CL (*> fc -L ) [ ) tpL+CL (y, l± ) 



dy d?lj 



d 2 J± Tt 



,QyCL 



2/(1-2/) lft? 7 dz i6^^L(^^)(^T)^ + CL(y,l^, 

(2.27) 



where ipL(x,k±) is the two-body wavefunction obtained by including only the ladder kernel V L 
while i/jl+cl( x , k±) and ipL(x,k±) are the wavefunctions including the cross-ladder kernel V CL 
and the self-energy correction term f(x,k±), respectively, in addition to V L . Likewise, we define 



see Appendix B for the full expression of Plowi Plowi Pmgh, an d Pffigh- 
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P£™ L+SB , P £1 CL+SB > p m£ L+SB i and P m g h L+SB b y addin S th e contribution from V s B . By 
adding the contribution from V HF to this, we may also define P£^ CL+SB+HF , P£^ CL+SB+HF , 
P i^h L+SB+HF , and P*j+° l +sb+hf _ Hcr ^ the factQr 167rTOlTO2 in the re l a tion between a and g 2 
has been compensated by the same factor included in the definition of the expectation values given 
by Eqs. (|2.11[) - (|2.15|) . We note that the spectrum of the bound-state is intimately correlated with 
the corresponding bound-state wavefunction. Thus, it is interesting to compute the probabilities 
Plow, -Plow, -Pffigh, and Pm g h with the specific contributions from the kernel and/or the self-energy 
correction and discuss the correlation of these wavefunction-related observables with the results of 
the spectrum. We present the probability ratios such as P£ gh /P£ w , P^ L /P^ L , Pm s J P L 
etc. in the following section of Numerical Results. 



Low ' 



3 Numerical Results 
3.1 Spectrum 

We present the results for the binding energy B = mi + m 2 — M as a function of the coupling 
constant a = g 2 / (MSirmim?) ■ Since we consider in general the m\ ^ ra 2 and A 7^ cases as well as 
the special cases of mi = 7712 or A = 0, we present the numerical results of representative examples 
(A/m 2 = 0,0.15,0.5, 1.0,2.5) with (3 = m 1 /m 2 = 1 and 5.446x 10" 4 in units where m 2 = 1. Whether 
the range of the interaction between the two constituents of the bound-state is short or long depends 
on the value of A: e.g. the Yukawa interaction with A =/= or the Coulomb interaction with A = 0. 
Also, the value of (3 is linked to the modeling of the bound-state: e.g. positronium or deutron system 
with [3 = 1, and hydrogen atom with f3 — 5.446 x 10~ 4 . 

Fig. [5] shows the case of (3 = 1, A = 0. The numerical results of the ladder, L + CL, L + CL + SB, 
and L + CL + SB + HF with and without the self-energy corrections obtained by the variational 
principle are presented and compared with the available numerical or analytic results. The non- 
relativistic result given by Eq. (|2.20j) is shown as the Balmer Line (thick grey solid line) in the far 
left of this figure, while the analytic result (Analytic L+self) given by Eq. (|2.19|) is shown as a thin 
grey solid line at the bottom of the figure. In the middle between these two curves (Balmer Line 
and Analytic L+self), our ladder result (thin black dot-dashed line) is compared with the previous 
numerical results (solid circles) provided by M. Mangin-Brinet and J. Carbonell [2^. Although our 
method of computation is different from theirs [33] , our ladder results are in good agreement with 
their ladder results. It lends confidence to our variational calculation. Just next to the ladder result 
on the right, the analytic result without the self-energy correction, i.e. without the last term in 
Eq. (12.191) . is shown as a thin grey dashed line. The analytic result agrees with the numerical result 
in the weak binding limit. Similarly, the analytic result including the self-energy effect (Analytic 
L+self) shown at the bottom of the figure is comparable to the corresponding numerical results 
(L+self) drawn as a thick black dot-dashed line just above the analytic result. Between the Balmer 
line and the ladder result, one can see the three curves of L + CL, L + CL+SB and L + C L+SB+HF 
from right to left. Also, between the ladder result and the L+self result at the bottom, one can 
see the three corresponding curves including the self-energy effect. These results indicate that the 
effect of the self-energy term is highly repulsive while the contribution from the cross-ladder shows 
a large attractive effect. The difference between the L+self and L + CL+self is more significant 
than the difference between the ladder and L + CL, which indicates that the attractive effect from 
the cross-ladder kernel is more significant when the self-energy corrections are included. We note 
that the particle-antiparticle creation/annihilation contribution from V HF as well as the stretched- 
box contribution from V SB yield attractive effects and become larger as the coupling constant a 
increases. Although the effect of the particle-antiparticle creation/annihilation is very small in this 
case (/3 = 1, A = 0), it appears more noticeable when the self-energy corrections are included as one 
can see in the difference between the thick black solid and thick black dashed lines. 

In Fig. [(J we consider the case of /3 = 1, A = 2.5 where we find that the particle-antiparticle 
creation/annihilation contribution from V HF becomes quite substantial. The contribution of V HF 
is much greater than that of V s B , with or without the self-energy corrections. In particular, with 
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the self-energy corrections, the V contribution becomes large enough to make the curve L + 
CL + SB + HF+self go above the ladder curve. We also note that the bound states can be formed 
only with a sufficiently large coupling since it is difficult to exchange a large mass particle with 
insufficient coupling strength. Moreover, due to the short range interaction in this case, the effects 
of the self-energy and cross-ladder seem to appear quite different when compared to those in the 
case of A = 0. Although we did not include the whole curve of L+self we note that a reaches 57.4 
when the binding energy B becomes 2m. 

In Fig. [3 the ladder (thin black dot-dashed), L + CL (thin black dotted), and L + CL + SB (thin 
black dashed) for the case of /3 = 1, A = 0.5 are compared with the corresponding numerical data 
(solid circles, squares, and stars, respectively) given by J. Carbonell and V.A. Karmanov 0, 
Our numerical results agree well with their results. This indicates that our variational wavefunction 
in Eq. (I2.17|) . based on the solution for the Hulthen potential, is reasonable for A ^ case. The thin 
black solid line for L + CL + SB + HF is also displayed in this figure to show the V HF effect. 

Fig. [5] presents the unequal mass case of (3 — 5.446 x 10~ 4 , A = 0. The ladder (thin black 
dot-dashed), L + CL (thin black dotted), and L + CL + SB (thin black dashed), and L + CL + 
SB + HF (thin black solid line) are shown in this figure. Here, we see that the particle- antiparticle 
creation/annihilation contribution from V HF is non-negligible although the exchange particle mass 
is zero (A = 0). All the results including the self-energy effects, such as L+self, L + CL+self, 
L + CL + SB+self, and L + CL + SB + HF+self, are not shown in the figure because they are all 
strongly suppressed in this case. 

Fig. [9] shows the case of f3 = 5.446 x 10~ 4 , A = 2.5. We find that the particle-antiparticle 
creation/annihilation contribution from V HF is enhanced most significantly in this case. We also 
note that the attractive effect of the cross-ladder contribution becomes significant when the self- 
energy is included as one can see from the L+CL+self L+CL+SB+self and L+CL+SB+HF+self 
curves. However, the L+self result is away from the other results and we note that a reaches 25.45 
as the binding energy B approaches to m. 

These results confirm that the self-energy effect is repulsive while the effect of the cross-ladder 
contribution is attractive. Comparing the results with and without the self-energy effect, we no- 
tice that the attractive effect from the cross-ladder contribution is more significant when the self- 
energy is included. Although the effects of the stretched-box and the particle-antiparticle cre- 
ation/annihilation are also attractive, their contributions are not as significant as the cross-ladder 
contribution. In particular, the effect of the stretched-box is further reduced when (3 ^ 1 and A ^ 0. 

Although the particle-antiparticle creation/annihilation contribution is not so significant when 
(3 = 1 and A = 0, it becomes non-negligible as (3 ^ 1 and/or A ^ 0. We find that the V HF effect 
is the most significant in the case of /3 ^ 1 , A ^ and also non- negligible for the (3 = 1, A ^ and 
[3 ^ 1, A = cases. 



3.2 Zero Binding Limit 

It has been indicated in [24, 3^, 35, 36 1 that for the non-zero mass of the exchange particle \, the 
relativistic ladder results in the weak binding limit B — > do not coincide with the non-relativistic 
result. We confirmed this result as shown in Fig. [101 Our results of ladder (thin black dot-dashed), 
L + CL + SB (thin black dashed), and L + CL + SB + HF (thin black solid line) are depicted 
for A = 0.15, 0.5, and 1.0, and compared with the available ladder (solid circles) and L + CL + SB 
(stars) results from Ref . [33] . Also our results of the Schrddinger equation with the Yukawa potential 
obtained by the variational principle (thick grey solid line) are compared with the corresponding 
results (diamonds) from Ref. [33[ . Our results are in close agreement with all of the available 
results from Ref. [33j | . As one can see in Fig. [101 the discrepancy between the non-relativistic 
limit and the zero binding limit becomes larger as the exchange particle mass increases. However, 
the cross-ladder and stretched-box reduce the discrepancy. We notice that the particle-antiparticle 
creation/annihilation contribution from V HF further reduces the discrepancy although it is not as 
significant as the cross-ladder and stretched-box. To show the more noticeable V HF effect, we 
present the A = 2.5 case in Fig. [TT] One may expect that the discrepancy would not be completely 
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Figure 5: /3 = 1, A = 0. The numerical results of the ladder, L + CL, L + CL + SB, and L + CL + 
SB + HF with and without the self-energy corrections obtained by the variational principle. The 
ladder and L+self are compared with the corresponding analytic results from Eq. {2T9}. The Balmer 
line is shown as the non-relativistic result from Eq. 112.201 1. The ladder is also compared with the 
previous data from M. Mangin-Brinet et al. 1 (M. Mangin-Brinet and J. Carbonell 33].) 




a 



Figure 6: /3 = 1, A = 2.5. The numerical 
results of the ladder, L + CL, L + CL + SB, 
and L + CL + SB + HF with and without 
the self-energy corrections obtained by the 
variational principle. 




a 



Figure 7: /3 = 1, A = 0.5. The ladder 
is compared with the corresponding ana- 
lytic results from Eq. l|2.21|l . The ladder, 
L + CL, and L + CL + SB are compared 
with the previous data. x (Data from J. 
Carbonell and V.A. Karmanov [24l | and 
personal communication with J. Carbonell 



11 



B 



1.0 


- - Ladder i 
L+CL i 

- - L+CL+SB / i 

- L+CL+SB+HF / i 




0.6 




0.4 




0.2 




0.0 





a 



Figure 8: @ = 5.446 X 1CT 4 , A = 0. The 
numerical results of the ladder, L+CL, L+ 
CL + SB, and L + CL + SB + HF without 
the self-energy corrections obtained by the 
variational principle. 
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Figure 9: j3 = 5.446 X 1CT 4 , A = 2.5. The 
numerical results of the ladder, L + CL, 
L + CL + SB, and L + CL + SB + HF with 
and without the self-energy corrections ob- 
tained by the variational principle. 



removed unless all the irreducible kernels are included in the relativistic bound-state equation. 

3.3 Wavefunction Renormalization and Probabilities of Lowest and Higher 
Fock-States 

In Subsection 2.4, we defined the probabilities of finding the lowest two-body Fock-state and the 
higher Fock-states with and without the self-energy correction. In particular, including the self- 
energy correction, we have renormalized the two-body wavefunction and defined the corresponding 
probability with the renormalized wavefunction (see Eq. (|2.25jl ). The probability to find the higher 
Fock-states with the self-energy effect was also defined by Eq. ()2.26|) . 

In this subsection, we present the numerical results of the probability ratios of Pmgh/-PLow without 
the self-energy correction and -Pmgh/^Low with the self-energy correction. In Figs. [T3] and [T31 we show 
the probability ratios for ^High/^Low ( thin black dot-dashed), P^ L / P lJw CL ( thin black dotted), 
Pm^ L+SB / P H CL+SB ( tkin black dashed), and P^ L+SB+HF /P£+ W CL+SB+HF (thin black solid 
line) as well as the probability ratios for P^( igh / 1 P[ ow (thick black dot-dashed), P^£ L /P^ L ( tbick 

black dotted), P^ L+SB /P£+ CL+SB (thick black dashed), and P^ L+SB+HF '/ ' P£+ CL+SB+HF 
(thick black solid line), for A = and 2.5, respectively, with (3 = 1. 

Wether or not the self-energy corrections are included, we note that the probability ratios become 
larger as the coupling constant a increases. This indicates that the sector of the higher Fock-states 
becomes more significant as a increases. The probability ratio of the ladder with the self-energy 
corrections is constantly lower than that of the ladder without the self-energy corrections due to the 
highly repulsive effect from the self-energy term. 

In the case of /3 = 1 and A = 0, as shown in Fig. [T21 the attractive effect from the cross-ladder 
term in P£+c L /P£+ w CL becomes so large that P F ^ L /P^ CL even crosses over P££ L /Pg? L 
around a = 1.5. This indicates that the attractive effect from the cross-ladder kernel is more 
enhanced when the self-energy corrections are included. This behavior of the enhanced attractive 
effect from the cross-ladder contribution with the inclusion of the self-energy correction is consistent 
with what we have observed in our spectrum calculation (see e.g. Fig. O and the corresponding text 
in Subsection 3.1). This consistency between the results of the wavefunction related observables 
and the mass spectra manifests the correlation between the wavefunction and the corresponding 
mass spectrum as the eigenfunction and the corresponding eigenvalue must be correlated in the 
bound-state calculation. 

Due to the highly repulsive effect from the self-energy, the probability ratio Pm g h/ P F ow (with self- 
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Figure 10: Zero binding energy region of the ladder, LD+CL+SB, and LD+CL+SB+HF results 
compared with the non-relativistic ones (Schrodingcr equation with Yukawa potential), ft = 1, A = 
0.15,0.5, and 1.0. 1 (M. Mangin-Brinet and J. Carbonell [jjjj.) 
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Figure 11: Zero binding energy region of the ladder, LD+CL+SB and LD+CL+SB+HF results 
compared to the non-relativistic ones (Schrodinger equation with Yukawa potential), ft = 1, A = 2.5. 
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Figure 12: /3 = 1, A = 0. The numerical results of the probability ratios of PHigh/^Low (without self- 
energy corrections) and fHigh/^Low (with self-energy corrections) for the ladder, L+CL, L+CL+SB, 
and L + CL + SB + HF. 



energy) is lower than -Pjiigh/ -^Low (without self-energy) indicating that the probability of finding the 
higher Fock-states is lower when the self-energy effect is included without any significant attractive 
channel such as the cross-ladder contribution. When the cross-ladder contribution is included, 
however, there is a dramatic increase in the probability of finding the higher Fock-states. Without 
the self-energy effect, the ratio P^ L /P^ CL is much higher than the ratio Pm gh /Pl ow - with thc 
self-energy effect, the increment of the ratio P^+£ ' L / ' ' L is more dramatic than that of the ratio 

^Hi^h/-^Low as tne coupling constant a becomes larger, even crossing over P££P V^Low^ around 
a =1.5. 

In the case of f3 = 1 and A = 2.5, as shown in Fig. [13l we see again that the effect of the 
cross-ladder contribution is highly attractive and enhances the probability of finding the higher 
Fock-states wether the self-energy corrections are included or not. However, the increment of the 
ratio ^Higjf V^Low Ci is not as dramatic as in the case of A = and the ratio P^ L /P^ CL does 
not cross over the ratio P^^ L / P^J^ L ■ This appears to be consistent with the reasoning that the 
probability of finding the higher Fock-states, including the very massive exchange particle of A = 2.5, 
is much less than the corresponding probability involving the massless exchange particle of A = 0. 

4 Conclusion 

We have developed a variational method for solving the scalar field model theoriy (mi ^ 7712, A ^ 0) 
up to second order in the coupling constant a. The solutions for the Schrodinger equation with the 
Hulthen potential can be taken as a reasonable variational wavefunction for the A ^ case. We go 
beyond the light-front ladder approximation and extend to the full LFD kernel including the ladder, 
cross-ladder, stretched-box, and particle-antiparticle creation/annihilation contribution (V ). We 
also obtain the light-front two-body equation including the self-energy corrections in the case of 
mi ^ 77?,2, A 0. Our numerical results of the light-front formalism agree well with the previous 
numerical data (2~3. |33|. They are also consistent with our analytic result in the weak binding limit. 
Our results between the coupling constant a and the binding energy B indicate that the self- 
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a 

Figure 13: ft = 1, A = 2.5. The numerical results of the probability ratios of Pfflgh /^Low (without self- 
energy corrections) and -PHigh/fLow (with self-energy corrections) for the ladder, L+CL, L+CL+SB, 
and L + CL + SB + HF. 



energy effect is repulsive while all other effects that we took into account, such as the ladder, cross- 
ladder, stretched-box, and particle-antiparticle creation/annihilation contributions, are attractive. 
Besides the ladder, the cross-ladder effect is most significant among the attractive contributions. 
The effect of the cross-ladder is more significant when the self-energy is included. 

Although the particle-antiparticle creation/annihilation contribution from V HF is not so signif- 
icant for j3 = 1 and A = 0, it becomes non-negligible as /? ^ 1 and/or A ^ 0. In particular, we find 
that the V HF effect is the most significant in the case of (5 ^ I , A ^ and also non- negligible in the 
/3 = l,A^Oor/3^1,A = case. In contrast to this behavior of the V HF contribution, the effect 
of the stretched-box V SB is more reduced as j3 ^ 1 and/or A ^ 0. 

We also note that the spectrum of the bound-state is intimately correlated with the correspond- 
ing bound-state wavefunction. We see that the attractive effect of the cross-ladder enhances the 
probability of finding the higher Fock-states whether the self-energy corrections are included or 
not. In the case of A = 0, without the self-energy effect, the ratio ^High V^Low ^ s mucn higher 
than the ratio ^High/^Low- Including the self-energy effect, the increment of the probability ratio 

^Higr? i /^Low C ' L i s morc dramatic than that of the ratio ^High/^Low as a i ncreases - This seems 
consistent with the result of the spectrum calculation exhibiting the more pronounced cross-ladder 
effect when the self-energy is included. As the exchange particle mass increases, however, the prob- 
ability of finding the higher Fock-states may be reduced as one can understand from the derivative 
of energy denominator involving the exchange particle mass. As such, in the case of A = 2.5, the 
increment of the ratio L / P^J^ ' L IS not as dramatic as in the case of A = 0. The probability of 

finding the higher Fock-states with the massive exchange particle of A = 2.5 is much less than the 
corresponding probability involving the massless exchange particle of A = 0. 
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Appendix A: Full Expression of the Kernel 

Corresponding to Figs. [U [3J and 01 respectively, the stretched-box (V s B ), cross-ladder (V CL ), and 
higher-Fock (V HF ) kernels are written in the form: 

V SB = ~, r £ Vi , (A.l) 

z(x - z)[z - y)(l - z) .f^ 2 

vCL = -, — w — 7n nE^> ( A - 2 ) 

z(x - z){z- y)(l - x - y + z) ^— J 

V HF = ^ — [Va + V B ] , (A.3) 

z[x — z)[z — y){l — x — y + z) 

where Vi, V 2 , V 3 , V 4 , V 5 , V 6 , V7, V 8 , Va, and V B are given by 

Vi ={M 2 - ^ +m i _ CJ-L -ki_) 2 + \ 2 _ jj + ml x 
x z — x 1 — z 

X (M 2 - ^ +m ' - (^-^) 2 + A2 _ (^-J±) 2 + A 2 _ ft+m| x 
x z — x y — z 1 — y 

x (m 2 - *L±m£ - ^-^ 2 + A2 - i±^l)-i , ( A.4) 

z y~z 1-y 

V 2 =(M 2 - S ± m i _ (k ± -j ± ) 2 + X 2 _ fcj +m| x 

z 2; — z 1 — a; 

x(M 2 - ^ +TO i _ - + A " _ 5 + m 2 )-i 
y z-y 1 - z 

x(M 2 ft + m| (fc±-J±) 2 + A 2 (J±-f±) 2 + A 2 fcj+m 2 ^ 
y a; — z z — y 1 — x ' 
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V 3 = - (M 2 - J± - "" 1 - ^ l± ~ J± > - *i +m 2 )- 



x(M 2 



r , /; - in] < ! -./_> t- * y fcj. - jj.) 2 + A 2 (j±-k ± -l±) 2 + mj x 

x — z 1 — x — y + z 

> i .u 2 - ^— ^ - > - ^ + ^ 2 )- 1 , (A.6) 



V A = — (M 
x(M 2 



y — z 1 — x — y + z 
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V A =(M 2 - ^ +m ' _ ( fc ^-.^) 2 + A2 _ (.U-fc±-^) 2 + mj ! 

y a; — z I — x — y + z 



x(M 



2 fcj +m\ _ jj+m 2 _ P ± +ml _ (j ± - k± - t ± ) 2 + m 2 x 
x z y 1 — x — y + z 



x(M 2 fcj_+m 2 (/±-J±) 2 + A 2 (ji-ii-li) 2 +m| ! 
a; y — z 1 — a; — y + z 

Vb =(Af 2 - S ± m i _ (fc-i ~ J-i-) 2 + A2 _ fcj + m l -i 
z x — z 1 — X 
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x (M 2 - 2+1^ - ^--^ 2 + A2 - i±^)-i . (A.13) 
z y - z 1 - y 

Appendix B: Full Expressions of Plow and Pmgh 

The full expressions of Pl ow and Pmgh without the self-energy can be given by the following 
integrations: 



dx erfcj 

x(l — x) 167T 3 



^L L ot CL - / I7^3^lfccL(^fci)| 2 , (B.2) 



P L -— = / ^^ + - + ^)| 2 , (B.3) 



i+CL= /_a(ay^+a 2 v^)\ 
High y aikf2 y 

2 / ^ ^J- / dy d 2 l ± )f -. , /WM - 

= -- 9 7 xTT^) 16^ J ^y)^ L+CL[x ' ^ \dNP)^ L+CL{vM 

4 fdxd 2 k ± f dy d 2 l ± f d 2 j ± t r.^Hi r 

High y aM 2 / 

2 f dx d 2 k ± f dy d 2 t ± , t , - ,fdV L \ , , 

= - 9 J xTT^yie^y ^^i^^+^ (x ' fci H^J^ +Ci+Sij(2/ ' ±} 

4 [ dxd 2 k ± f dy d 2 ^ f d 2 J ± f - 



(B.5) 



4 / dx d 2 k ± f dy d 2 l ± f d 2 j ± t - N ^dV SB 



(B.7) 
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P L+CL+SB+HF _ _ 
^High ~~ 1 



d(aV L + a 2 V CL + a 2 V SB + a 2 V HF ) 



dx d 2 k\ 



dM 2 
dy d 2 l± | 



.t(1-x)16^ 3 J y(l-y) iq^l+cl+sb+hf^,-^; y m 



( dV L ' 

(^-1-) ( flI72 ) ^L+CL + SB+HF{y,ll_) 



dx d 2 k± 



f dy d 2 ly_ f d 2 j 1 _ ;t , 

J y(l - y) 16^3 J dZ 16n ^L + CL + SB + HF^^^ 



x(l — x) 16-7T 3 J y(l — y) 167r 3 
dx d 2 k± f dy d 2 l± 



r , fdV CL 
fell 



V <9M 2 



1pL+CL + SB+HF(yJ±) 



d 2 ]±_ 

(1 y(l - y) 16tt 3 ./ ^16^^ ' ^ ' " ; " 



4 /" dx d 2 k± r 
~ 9 J x(l - x) 16tt 3 J 



QySB \ 

^V'l+CL+SB+HF^^-L) ( -g^r) 4>L+CL + SB+HF(y,l±) 



dx d 2 k i 



dy cPl± 



x(l — x) 16-7T 3 J y(l — y) 16-7T 3 



dz 16ir 3 ^l+cl+sb+h I 



(x, kj_) 



dV 



HF 



dM 2 



lpL+CL+SB+HF(y, l±) , 

(B.8) 



where ip(x, k±) is the two-body wavcfunction optimized by the variational parameter for the corre- 
sponding kernel and coupling constant a. Likewise, the full expressions of Plow and Pmgh with the 
self-energy can be given by the following integrations: 



pL 
Low 



pL+CL 
Low 



/dx d 2 k± 
X(l — x) 167T 3 



pL+CL+SB _ 
Low 



pL+CL + SB+HF _ 
Low 



dx d 2 k± 



x(l — x) 167T 3 
dx d 2 k 



_ am 1 m 2 df(x,k±) 

[ 7T ' dM 2 

amim 2 df(x, k±) 



l-( 



X{1 — x) 1671" 3 

dx d 2 kj_ 

x(l — x) 167T 3 



,r ; dM 2 

OCTTl\Tfi2 . Of(x, k±) 

5? ^ <9M 2 

amim2 . df(x, k±) 
5? ' dM 2 



\M*,k ± )\ 2 , 

\lJ>L+CL(x,k±)\ 2 ■ 



(B.9) 

(B.10) 
(B.11) 



IV'L+CL + SB+Hf (x, k ± )\ 2 , (B.12) 



p 



High 



= -.9 



<9M 2 

— x) 167T 3 



dy d 2 l 
y(l-y) IStt 3 ^ 



-> ( dV \ — 

-4>{(x, k ± ) [-q^J i>L(y, T±) , 



(B.13) 



pL+CL _ 
^High - 



-.9 



dM 2 



/" dx d 2 k± f 1 
J xTl - x) 16tt 3 J ~y~(i 



x(l — x) 167T 3 J y(l — y) I671" 3 
dx d 2 k± f dy d 2 l± f , d 2 j± - 



^l+cl( x ^-l) ( T^H ) i>L+CL(y,l±) 



dz- 



x(l - x) 16tt 3 J y(l - y) 16tt 3 J 16tt 3 ^ l + cl 



, fdV CL \ 
X x ' fc J-) -5775- *pL+cL(y, l±) , 



V <9M 2 



(B.14) 
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pL+CL+SB _ I _ 



0(aV L + a 2 V CL + a 2 V SB )\ 



High \ q M 2 I 

2 f dx d 2 k ± f dy d 2 l ± 7f -. /9^\- - 

= ~ 9 J i^T^lfe? 7 J^^^UL + Sfl(*,*±)^ > )^ i+ CL + SB(y,ix) 

4 f dxd 2 k ± f dy d 2 ^ f d 2 J ±7f j. /c^ Ci \ - - 

4 f dxd 2 k ± f dy d 2 l ± f rf 2 J ±7t r,/W SB \j , 

(B.15) 

Hi « h \ 9M 2 / 

= ~ g J xJl=xjT6^ J ^V)16^^l + sb + hA^) ^—)^ L+CL+SB+HF (y,l ± ) 

4 fdxd 2 ^ fdyd 2 ^ f d 2 j ±Ji t (8V CL \ ~ r 

4 f _Jxd 2 k ± f dy d 2 l ± f d 2 j ± ~i t fdV SB \ ~ f 

4 [ dxd 2 k^ fdyd 2 ^ f d 2 J ±7t -. ( dV HF \ j - 
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